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The flavor-conserving non-leptonic weak interaction can be studied experimentally through the 
observation of parity violation in nuclear and few-body systems. At hadronic scales, matrix elements 
of parity- violating four-quark operators ultimately give rise to the parity violating couplings between 
hadrons, and such matrix elements can be calculated non-perturbatively using lattice QCD. In this 
work, we investigate the running of isovector parity-violating operators from the weak scale down 
to hadronic scales using the renormalization group. We work at next-to-leading order in the QCD 
coupling, and include both neutral-current and charged-current interactions. At this order, results 
are renormalization scheme dependent, and we utilize 't Hooft-Veltman dimensional regularization. 
The evolution of Wilson coefficients at leading and next-to-leading order is compared. Next-to- 
leading order effects are shown to be non-negligible at hadronic scales. 

PACS numbers: 12.38.Bx,12.15.Mm,12.38.Cy 



I. INTRODUCTION 

The tremendous progress made in understanding weak 
interactions led to the development of the Standard 
Model of particle physics. Many corners of the Standard 
Model have been tested to very high precision through ex- 
periment. The hadronic neutral weak interaction, by con- 
trast, is the least constrained of all Standard Model cur- 
rents. Such an interaction is difficult to probe; because, 
even while largely shielded from electromagnetic contri- 
butions, it is easily overshadowed by dominant strong 
contributions. One must look to processes that vio- 
late strong interaction symmetries to isolate the hadronic 
weak interaction. Flavor-changing neutral currents are 
suppressed in the Standard Model due to the so-called 
GIM mechanism [l[. To see the hadronic weak interac- 
tion through flavor-conserving neutral currents, one must 
look to isovector hadronic parity violation^ Recent ex- 
perimental advances give reason for a concerted effort to 
study hadronic parity violation theoretically using QCD. 

Experimentally hadronic parity violation has been 
sought in nuclear reactions. In fact, the first search for 
hadronic parity violation in the nucleon-nucleon inter- 
action was published the same year (1957) that parity 
violation was discovered in the weak decay of nuclei [2|. 
Since then, hadronic parity violation has been observed 
in a variety of nuclear reactions, for example, the first 
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1 The charged-current interaction makes contributions to isovec- 
tor hadronic parity violation at the weak scale that are doubly 
Cabbibo suppressed, i.e. oc \V US | 2 /| V u d\ 2 ~ 0.05. As we work to 
next-to-leading order in the QCD coupling, we include such sub- 
dominant contributions in our analysis. At hadronic scales, we 
find that coefficients of strangeness containing |AJ| = 1 operators 
are modified by the partity violating charged-current interaction 
at the 10-20% level. 



conclusive evidence came from radiative neutron capture 
on 181 Ta and a substantial (~ 10%) parity- violating 
effect was observed in 139 La by using orders-of-magnitude 
enhancement produced from a near degeneracy of oppo- 
site parity states that mix under the weak interaction [J] . 
A more thorough summary of these and other experi- 
ments appears in the reviews [H, Q. Recently bounds 
on hadronic parity violation in few-body systems have 
been obtained by measuring the spin rotation of a neu- 
tron traveling through 4 He [7], and through the measure- 
ment a parity-violating spin asymmetry in the reaction 
n + p d + 7 @. 

Using the experimental data on hadronic parity viola- 
tion in nuclear reactions to probe the hadronic neutral- 
current of the Standard Model is an ambitious under- 
taking. Interpretation of the experiments is usually car- 
ried out using the meson-exchange model of Desplanques, 
Donoghue, and Holstein (DDH) In this model, vari- 
ous parity- violating meson-nucleon couplings give rise to 
the parity- violating nucleon-nucleon interaction. This in- 
teraction then generates hadronic parity violation in nu- 
clei. The parity-violating couplings are unknown model 
parameters (although estimates and broad natural ranges 
are provided by DDH). Assuming the DDH framework, 
constraints on the couplings can be extracted from ex- 
periments. While the results agree within the reason- 
able ranges, constraints from different parity-violating 
experiments are not consistent. In this framework, it 
is challenging to determine where discrepancies may lie: 
in model assumptions buried in nuclear structure calcu- 
lations, in model assumptions about the parity-violating 
nucleon-nucleon force, in dynamical effects due to the 
non-perturbative nature of QCD; or perhaps, although 
very unlikely, in the non-leptonic weak interaction itself. 

Recently great strides have been made to remove model 
assumptions about the parity-violating nucleon-nucleon 
force [l(J [H| • For few nucleon systems at very low en- 
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ergies, one can use effective field theory techniques to 
parametrize the most general parity-violating nucleon- 
nucleon interaction based solely on symmetries of the 
system^ In this framework, there are five unknown 
parity-violating couplings [13j that can in principle be 
determined from few-body experiments, see, for exam- 
ple, There is much anticipation that, in the 
not-too-distant future, values for some of these couplings 
will be constrained due to the remarkable experimental 
program dedicated to measuring hadronic parity viola- 
tion in few-body systems. It is not clear when complete 
information about the parity-violating nucleon-nucleon 
interaction will be known. Even if the parity-violating 
nucleon-nucleon force is fully mapped out, the connec- 
tion of its terms to the weak interaction, moreover, can 
only be done through non-perturbative QCD. 

In lieu of first principles information about matrix 
elements of these hadronic-scale, parity- violating op- 
erators, symmetry arguments have been used to ana- 
lyze sources of parity violation at the hadronic level. 
Most notably, there was a mo del- independent analysis of 
pion-nucleon couplings using chiral symmetry considera- 
tions [llj]. This analysis systematically enumerated the 
low-energy parity- violating couplings between pions and 
nucleons, as well as pions, photons and nucleons. Siz- 
able contributions in the |A7| = 1 sector (for example, 
to the isovector parity-violating pion-nucleon coupling 
^vNn) were argued to arise from operators involving an 
ss-pair. Using chiral perturbation theory, the one-loop 
chiral corrections to were later determined fliij . 

This particular coupling is believed to be phenomenolog- 
ically important, as it gives rise to the long-range part 
of the parity-violating nucleon-nucleon force. Theoreti- 
cians have pointed to ways this coupling might pos sibly 
be measured through high-statistics experiments [2CJ-(23| . 
Anticipating lattice calculations of the parity-violating 
pion-nucleon coupling, a chiral perturbation theory anal- 
ysis was presented using partially quenched chiral per- 
turbation theory [24] H 

Lattice QCD computations can determine hadronic 
parity violation from first principles. Quite recently the 
first lattice QCD study of nuclear parity violation ap- 
peared [26| . In this pioneering work, the isovector chan- 
nel was the focus, and a signal for the parity-violating 
pion-nucleon coupling, h^ N , was obtained. As with 
any lattice computation, there are a number of sys- 
tematic errors which must be controlled to make con- 
tact with phenomenology. As further refinements to the 



method are made, we can expect to see precision infor- 
mation about parity- violating hadronic couplings coming 
from lattice QCD computations]! Complete information 
about hadronic parity violation in few-nucleon systems 
may not come solely from experiment or theory. A com- 
bination of the two will thus be essential to provide input 
for hadronic parity violation in nuclear many-body cal- 
culations. To this end, we study how well the sources of 
parity violation are known at hadronic scales. A leading- 
order QCD analysis to this effect was presented some 
time ago in [29(. Here we focus on the isovector parity- 
violating operators, and include higher-order effects. 

Starting from the weak currents in the Standard 
Model, we determine the isovector parity-violating in- 
teraction at weak scales to next-to-leading order accu- 
racy in the strong coupling. We include contributions 
from both neutral-current and charge-current interac- 
tions. The resulting isovector parity-violating four-quark 
operators are then evolved down to hadronic scales us- 
ing the QCD renormalization group at next-to-leading 
order. To achieve this, we borrow heavily from the two- 
loop formalism extensively developed in the context of 
flavor-changing weak decays, see [301 ] for a thorough re- 
view. We find that next-to-leading order QCD effects on 
isovector parity- violating operators are non-negligible at 
hadronic scales. Renormalization of parity-violating op- 
erators, moreover, is a key ingredient that must be un- 
derstood in order to compute parity-violating hadronic 
observables on the lattice. Although we carry out this in- 
vestigation using a perturbative renormalization scheme 
defined in the continuum, our findings are nonetheless 
useful as they can serve as a gateway to convert to other 
schemes. 

Our presentation has the following organization. First 
in Sec. [Ill we enumerate our operator basis for isovector 
parity violation. This is carried out in three effective field 
theories: one valid just below the weak scale, another at 
intermediate scales, and finally the target effective field 
theory valid at hadronic scales. In Sec. IIII1 we deter- 
mine the isovector parity- violating operators at the weak 
scale to next-to-leading order in the strong coupling. Re- 
sults throughout are renormalization scheme dependent, 
and utilize dimensional regularization with the 't Hooft- 
Veltman scheme, in which 75 is non-anticommuting [Sir- 
|33j]. Next in Sec. lIVl the Wilson coefficients are run down 
to hadronic scales using the two-loop anomalous dimen- 
sions. Finally in Sec. |Vj we summarize our findings, and 
give an outlook to future work. 



2 The effective theory in question is referred to as the pion-less 
effective theory, see 1131 . This name obscures the fact that it is a 
universal low-energy theory valid for systems with large scatter- 
ing lengths, whether they be few-body nuclear systems without 
pions, atomic systems, etc. . 

3 The extension of electroweak charges of quarks advocated 
by [S^l , however, makes the computation difficult to use in prac- 
tice, see [25ll for further details. 



Lattice computations of hadronic parity violation are cur- 
rently at quite an early stage, and are considerably demand- 
ing. Decades worth of effort has been expended in an analogous 
pursuit, namely the study of matrix elements of flavor-changing 
weak interactions using lattice QCD. Quite recently significant 
progress has been made towards complete computations, for ex- 
ample, in the |A5| = 1 sector, see [27l, |2S( . We are optimistic 
that such success will carry over into the study of hadronic parity 
violation in lattice QCD. 



3 



II. OPERATOR BASIS 

Parity violation in nuclear and few nucleon systems 
arises from the non-leptonic weak interaction. In the 
isovector channel, the leading component stems from ex- 
change of a Z° boson, with a sub-leading ~ 10% contribu- 
tion from the exchange of bosons. At energy scales 
below the weak interaction, we can replace the weak bo- 
son exchanges by an effective theory containing contact 
interactions. 

The parity-violating, four-quark operators in the ef- 
fective theory undergo renormalization from gluon ra- 
diation. This radiation is important for evolving the 
operators from the weak scale down to hadronic scales, 
where their matrix elements can be computed using lat- 
tice QCD. As is standard, we must run the five-flavor 
theory down to the bottom quark mass, and match onto 
a four-flavor theory. This theory is then run down to 
the charm quark mass, where it is matched onto a three- 
flavor theory that is valid at hadronic scales. We first 
enumerate a closed basis of operators for each effective 
theory. 



A. Five-Flavor Effective Theory: m b < fj, < Mz 

At scales above the bottom quark mass, we have a 
five-flavor effective theory. The isovector parity- violating 
Lagrangian can be written in the following form 

(5) 

Notice that the effective Lagrangian C P y is independent 
of the renormalization scale /i. The evolution of the Wil- 
son coefficients Cf([i) is exactly compensated by the 

evolution of the operators O f (/x) . There are eight such 
operators. 

A complete basis of isovector parity-violating opera- 
tors in the five-flavor effective theory can be grouped into 
three sets. The first two sets, 

5 

of = {u lix u - d llx d) L X>7"<7)£ - {L R], 
i 

5 

of = [u llx u - d llx d] L Y^rfih ~{l^r}, 

q 

5 

of = (H 7/1 « - d^d) L YJjtfq)R ~{L^ R}, 
q 

5 

of = (uj^u - d lf >d] L J}^?)* ~ {L <* R}, (2) 

i 



and 

of = (u^u - d lf ,d) L {s^s - c^c + b^b) L 
-{Lh R}, 

of = (u^u - d^d] L [sj"s - c^c + b^b) L 

R}, (3) 

arise from Z° exchange. The subscripts L and R repre- 
sent the chiral components of the fermion bilinears, so 
that (4>'Y^)l : r = 4>l,r1h4>l,r, with ip L>R = V l ,r4>, 
and the chiral projectors are defined in the usual way, 
~Pl,r = 5(1 -F 75)- I n eacn case, the (tpj^t/j) notation in- 
dicates a color-singlet contraction, 4' a l^' i l }a - On the 
other hand, the mixed bracket notation denotes the fol- 
lowing color contraction, 

(M] [fiffi =Y^ a ^ h ^^ a - ( 4 ) 

a . b 

The chiral basis we employ is different than that chosen 
for isovector parity violation in [29] . For contributions 
arising from Z° exchange, the chiral basis enables us to 
eliminate two unnecessary operators in each effective the- 
ory. This point will be discussed further in Sec. IIIII 

Exchange of PF-bosons makes a smaller ~ 10% con- 
tribution to isovector hadronic parity violation, and was 
not considered in (29j . This contribution is included here 
(see Sec. MI Bl below). and the operators arising from W 
exchange are chosen to be 

Of = {uj^u - d^d) L (s^s - &fc) L - {L o R}, 
Of = (u^u - d^d] L [s^s - Srfc) L - {L o R}. 

(5) 

Taken together, the first two sets of isovector parity vio- 
lating operators are closed under QCD renormalization. 
Additionally the third set is also closed under QCD renor- 
malization. 



B. Four-Flavor Effective Theory: m c < (x < rrit 

Integrating out the bottom quark, we have a four-flavor 
effective theory. The isovector parity-violating effective 
Lagrangian is written similarly to that above, namely 



4^ = ^5>|%)^ 4 V)- (6) 
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There are now six operators in the effective theory; and, 
we choose the basis 

4 

of = (n^u ~ d^d) L X^gU - {L R}, 

1 

4 

of = (u^u - ^liEl^^ ~ ^ i?}, 

9 
4 

i 

4 

of = (u^u - d^d^YslQl^M -{L^ R}, (7) 

and 

Of = (u^u - d lfi d) L (sj»s - c^c) L - {L <-> 
C^ 4) = (u 7m u - d^d] L [s^s - Srfe) L - {L o i?}. 

(8) 

At these scales and consequently below, the contribu- 
tions from neutral-current and charged-current interac- 
tions become indistinguishable (up to sub-percent con- 
tributions we neglect, see Sec. IIIIBI below). Further- 
more, the operators shown in Eq. ([7]). and those oper- 
ators shown in Eq. (|SJ) each form a closed set under QCD 
renormalization. 

C. Three-Flavor Effective Theory: Aqcd < /U < m c 



and 

Of = (uj^u - d^d) L (s^s) L -{L^ R}, 

Of = {u^u-d^d] L [s^s) L -{L O R}. (11) 

These sets of operators mix with one another under QCD 
renormalization. Furthermore, not all of the operators 
are independent, because a Fierz transformation yields 
the relation 

Of = Of -Of + Of. (12) 

We will treat this operator as independent, however, and 
use the Fierz constraint as a consistency check for the 
anomalous dimension matrices. This is only possible be- 
cause the 't Hooft-Veltman scheme respects Fierz trans- 
formations (34T |. 



III. ISOVECTOR PARITY VIOLATION AT THE 
WEAK SCALE 

Having enumerated an operator basis for isovector 
hadronic parity violation, we now determine the parity- 
violating component of the non-leptonic weak interaction 
at the scale of weak interactions. We consider tree-level 
exchanges of weak vector-bosons, as well as the one- loop 
QCD radiation required at next-to-leading order. 



Finally at scales below the charm quark mass, we 
match onto a three-flavor effective theory. This is the the- 
ory used in lattice QCD simulations, where the hadronic 
scale roughly corresponds to the inverse lattice spacing 
fi ~ a -1 . The lattice regularization is, of course, a differ- 
ent scheme than the one considered here. The three- 
flavor isovector parity-violating effective Lagrangian is 
written the same way as above 

C% = ^£,Cy)(M)0?>to. (9) 

For the six terms of this three-flavor effective theory, we 
choose the operator basis 

3 

Of = {u lix u - d llx d) L X>7"<7)i - {L R], 



O 



(3) 



(u-y^u- d-f f j,d) L ^2[qj t "q)L - {in R}, 



O 



(3) _ f„. 



{u^u - d ltl d) L Y^^^R ~{L^ R}, 



of = (u^u - d^d} L J2[qri)R -{L^ R}, (10) 



A. Isovector Parity- Violating Neutral Current 

1. Tree Level 

In the Standard Model, the hadronic neutral current 
has the form 



J,f = 



1 



M cos 6 W 



sm 



lw (tpliiQip)] 



where -0 is the isodoublet of quark fields 



(13) 



, with U = c , and D = \s . (14) 



t 



The third component of weak-isospin, T 3 , is given by 
T 3 = | diag(l, — 1), and the electric charge matrix has 
the form Q = ± +T 3 . 

At the scale /i = Mz and below, we can integrate out 
contributions from the top quark to arrive at a five-flavor 
effective theory. In this theory, furthermore, the Z° can 
be simultaneously integrated out, leaving only local four- 
quark interactions. At tree level, this procedure results in 
the | AI | = 1 parity- violating effective Lagrangian given 
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in [29|, namely 



2. One-Loop Corrections 



-(5) 



Gp 
V2 



1 - 5 

- sin 2 8 w (uu - dd) A ^(qq)^ 



+ ( — — sin 9w ) (m« — dd)A(ss — cc + 66)y 



+ ( — — sin 6w ) (uu — dd)v(ss — cc + 66)^ 



(15) 



Above, denotes the Fermi coupling constant, Gp — 
V / 2g 2 /8M^, and the subscripts V and A represent the 
vector and axial- vector quark bilinears, ip^f^ip and ipl^lb, 
respectively. Notice we use a subscript on to denote 
that both the coefficients and operators are taken at the 
scale \i = Mz- 

The operators appearing in C^L Mz are not closed un- 
der QCD renormalization, and at scales n < Mz other 
operators will be radiatively generated. A complete set 
of such operators has been enumerated in Sec. UU One 
should note that the non-penguin operators in Eq. (|15[) 
come in exactly the combination A ® V + V ® A. In the 
chiral basis, this combination is proportional to the op- 
erator structure Vl <S) Vl — Vr ® Vr, where Vl,r refers to 
the fermion bilinear (ip^f^^L^- The orthogonal combi- 
nation of operators has the structure Vl®Vr — Vr ® Vl ■ 
Non-penguin operators of the latter type are not gener- 
ated from QCD radiative corrections at two-loop order, 
which is borne in by the block-diagonal structure of Taq 
in Eq. (j57| below. This feature will persist to all orders in 
regularization schemes that preserve chirality 35]. The 
operator basis we employ thus contains the minimal set 
of operators necessary to describe isovector parity vio- 
lation. Compared to the basis chosen by [20|, we have 
two fewer operators from Z° exchange in each effective 
theory. Taking into account the Fierz constraint, there 
are only five operators at hadronic scales. 

While using an operator basis with unnecessary oper- 
ators has no physical consequences, the computation of 
additional matrix elements using lattice QCD can con- 
sequently be unnecessarily costly. There is no computa- 
tional advantage to using A®V+V®A operators instead 
of each one individually, however, as we know of no alge- 
braic way to combine the two sets of required quark con- 
tractions. There is an advantage to using V<S>A and A®V 
operators over Vl®Vl— Vr®Vr and Vl®Vr—Vr®Vl op- 
erators; because, in the latter case, parity violation arises 
from the difference of parity-violating hadronic matrix 
elements, which is presumably statistically noisier than 
in the former case. Our minimal operator basis is thus 
economical only for perturbative calculations. We will 
present our final results for Wilson coefficients in the op- 
erator basis of [IH, which is practical for lattice compu- 
tations. 



As we work at next-to-leading order, we must addition- 
ally consider the one-loop QCD corrections arising from 
integrating out the top quark and Z° boson. At this or- 
der, the Wilson coefficients are determined by comparing 
one-loop graphs in the full and effective theories. 

In the full theory, the next-to-leading order isovector 
parity- violating Lagrangian can be written as 



-(5) _ 



r< 6 

(_T R \ ~\ 



1 = 1 



A i {M z )0f ) {Mz) 



(16) 



where the amplitudes Ai(Mz) are computed up to one- 
loop order, and are thus of the form 



A i (M z ) = (c ) i + ^^-(c 1 ) i 



■in 



(17) 



The (co)i are determined at tree level, and can be read 
off from Eq. (fT5|) . namely 



(co)i = Ma = -3 sin 2 6 W = -0.0771, 

(c ) 5 = -1 + 2 sin 2 9 W = -0.5376, (18) 



with all others zero. The (ci)j arise from the one-loop 
corrections, as we detail below. 

At the scale /z = Mz, the effective theory must re- 
produce Eq. (fT6|) at one-loop order. This is achieved 
by specifying values for the Wilson coefficients at this 
scale. At tree level, the Wilson coefficients are obviously 
given by C^ 5 \Mz) — (cn)i- At one-loop order, however, 
these coefficients are renormalizcd. After regularization 
and subtraction, the finite renormalization in the effec- 
tive theory will have the forrr@ 

C =^t M (% + 0f(Mz). (19) 

To reproduce the full theory, we must have 



cf\Mz) = (c ) i + 



a s (M z ) 



4tt 



(ci)i - V(c ) 



3 



(20) 

To compute the Wilson coefficients at ncxt-to-lcading 
order, we must then determine the finite terms (ci), in 
the full theory, and the matrix r in the effective the- 
ory. The former arises from the one-loop QCD correc- 
tions to Z° exchange, shown in Fig. [TJ whereas the latter 



Technically to derive Eq. H19H , the matrix r is computed by talk- 
ing certain matrix elements of the operators, (CJ\ (/x) }. One 
must then isolate terms that arise from the operator renormal- 
ization from those that arise from the particular external states 
chosen to compute the matrix elements. 
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FIG. 1. One-loop diagrams in the QCD renormalization 
of Z° exchange. The straight lines are quarks, while the 
curly lines are gluons, and dashed lines represent Z° bosons. 
The wave-function renormalization contributes but is not de- 
picted. 




FIG. 2. One-loop diagrams in the QCD renormalization 
of parity violating four-quark operators. Diagrams (a')— 
(/') show the renormalization of current-current interactions, 
while diagram (g') is a QCD penguin. Filled circles denote 
the four-quark operator, while all other diagram elements are 
as in Fig. [1] The wave-function renormalization contributes 
but is not depicted. 



arises from the one-loop corrections to the operators in 
Eq. (fT5|) . shown in Fig. [2] The finite contributions to 
these diagrams are renormalization scheme dependent. 

To perform these computations, we use dimensional 
rcgularization with MS subtraction in the 't Hooft- 
Veltman scheme for 75. Accordingly there are a total 
of d = 4 — 2e matrices 7^, that satisfy {7^,7^} = 2g^, 
where g^ v is the d-dimensional metric tensor. The defi- 
nition of the matrix 75 that appears in four-dimensions 
is retained in d-dimensions, so that 75 = — «7o7i7273- 
Thus if we write the Dirac matrices in terms of four- and 
(— 2e)-components, 7^ = 7 M + then 75 anti-commutcs 
in four-dimensions, {7^,75} = 0, and commutes in the 
remaining dimensions, [7^,75] =0. In extending the 
definition of chiral bilinears of the form ^iji) l ,r to 
d-dimensions, we must use the symmetric projection 

(lp1p1p)L,R = ^Vr^lJi/Pl^ = ^%Vl,riP- (21) 

Having spelled out almost all the details of our renor- 
malization scheme, we now proceed to match the full 
theory onto the effective theory at one- loop order in ot s . 
In performing this matching, it is trivial to integrate out 
the top quark, as there are no diagrams with an internal 
top quark appearing at this order. There is a top quark 
loop in the Z° self-energy, however, this contribution is 
automatically taken into account by using the physical 
value for the mass Mz- 

Ultraviolet divergences present in the one-loop QCD 
corrections to the current vertices [diagrams (a) and (6) 
of Fig. Q] are exactly cancelled by the wave-function 
renormalization. Ordinarily finite contributions are ad- 
ditionally cancelled; however, in the MS subtraction 
scheme with 't Hooft-Veltman rcgularization, the axial- 
current vertex receives a finite renormalization [36| . 
Specifically the renormalized axial-vector current has the 
form 

A„ = (l + ^CV) @ l %T 6 'h) > (22) 

with Cf — h(Nc — l)/-^cj and for any two quarks, ipi 
and tp2- Such renormalization implies a non- vanishing 
two-loop anomalous dimension for the axial-vector cur- 
rent in this scheme. Because the operators in Eq. ((15]) 
are each products of color-singlet bilinears, a similar fate 
befalls the analogous vertex corrections in the effective 
theory [shown as diagrams (a') and (b') in Fig. [2]. The 
ultraviolet divergences in the effective theory are can- 
celled by wave-function renormalization, but finite con- 
tributions are left. The difference of these finite con- 
tributions between the full and effective theories van- 
ishes. In this scheme, there is no 0(a s ) contribution in 
Eq. (|20p from matching the one-loop current vertex dia- 
grams. This scheme, however, has undesirable features, 
such as the mixing of left- and right-handed components 
due to Eq. ((22]) . and the appearance of the axial- vector 
current anomalous dimension in the evolution of the Wil- 
son coefficients. To preserve chirality, we augment the 
't Hooft-Veltman scheme by imposing a multiplicative 
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renormalization on the axial-vector current. The multi- 
plicative factor is chosen to force the renormalized axial- 
vector current to have the form 



A„ 



(V>i7m75^ 2 ) + 0(a 2 s ). 



(23) 



This augmented scheme is commonly employed to pre- 
serve chirality, see, e.g., [37|, HI]- To impose the aug- 
mented scheme, we recompute the current- vertex dia- 
grams in the full theory with the additional multiplicative 
renormalization. The next-to-leading order contributions 
from diagrams (a) and (6) in Fig. [T] accordingly vanish. 
Diagrams (a') and (&') in the effective theory now give 
non-vanishing contributions to the Wilson coefficients, 
upon taking the difference in Eq. (|20|) . 

In matching the penguin contributions in the full and 
effective theories, diagrams (g) and (g') of Figs. Q] and 

[3J respectively, we note that only the operator has 
a non- vanishing penguin contraction in the effective the- 
ory. The vanishing of penguin contractions for the other 
operators can be easily demonstrated. For the opera- 

(5) — 

tor 5 , the color-singlet bb cannot mix with the gluon; 

(5) 

and, for the operator Og , the left- and right-handed 
field cannot be contracted at one-loop order [an all or- 
ders argument is given after Eq. (fTUj) below]. Carefully 
evaluating the penguin diagrams in the full and effective 
theory, we find exactly the same contractions are involved 
as the penguins considered in [33] ■ In 't Hooft-Veltman 
regularization, moreover, the two diagrams are described 
by the same function, and their difference exactly cancels 
when the scale /i = Mz is specified in the effective the- 
ory. The remaining next-to-leading order contributions 
to Eq. ([20]) arise from comparing diagrams (c)-(f) in the 
full theory, to (c')-(f') in the effective theory. The for- 
mer diagrams are finite, while the latter require regular- 
ization and subtraction. Identifying the renormalization 
scale of the latter diagrams with Mz, moreover, ensures 
that there are no logarithms in the matching conditions, 
Eq. QQJ. 

Complete one-loop results for diagrams (a)-(f) in the 
full theory are described by (c\)i, whose non-zero entries 
are given by 

(Cl)l = ~(ci)2 = +^(c )l, 
(Cl)3 = "^( C l)4 = -^( c 0)3, 

(ci)s = -^(ci)e = +^(c )5- 



(24) 



This result includes the multiplicative renormalization 
of the axial- vector current. In the effective theory, the 
finite contributions of diagrams (a')-(f') in Fig. [2] are 
described by the matrix r in Eq. (|19[l . The relevant rows 
of r required for matching are 



r lj = 


(7 


-5 















r 3j = 


(0 





3 


7 









r 5j = 


(0 











7 


~ 5 ), • 


(25) 



Combining the results for (ci)j and r, we can match 
the two theories at fi = M z , and deduce the Wilson co- 
efficients, cf ] {M z ). We find 



C^{Mz) 



C { i\M z ) 

cf\Mz) 
cf\M z ) 



C^(Mz) 



1 - 



13 a s (M z ) 



2 4tt 
7 a s (M z ) 
2 4tt 



(co)i 



7 a s (M z ) 
2 4tt 



(co)i 



(co)i 



U a.(M 2 
2 4tt 
. 13a s (M z ) 



(co)i 



2 4tt 
7 a s {Mz) 
2 4tt 



(co)s 



(co)e 



(26) 



The next-to-leading order corrections give rise to a ~ 5% 
contribution to the input values of the Wilson coeffi- 
cients. 



B. Isovector Parity- Violating Charged Current 

The charged-current interaction also makes contribu- 
tions to isovector hadronic parity violation. The non- 
leptonic charged current has the form 



J, 



1 

71 



U^VD) 



(27) 



with = 



(Jj )t. The unitary matrix V is the CKM 
matrix (which should not be confused with the subscript 
V denoting the vector spinor contraction). At scales be- 
low Mw, we can integrate out the top quark and W- 
bosons simultaneously. At tree level, the resulting five- 
flavor isovector parity-violating effective Lagrangian has 
the form 



,/ (5) 
"u=M« 



2Gf 



(U^VD) L (pfV*U) L -{L^R} 

(28) 



Here the top quark has been removed from the spinor U, 
and the subscript on C' ^ is used to denote that both the 
coefficients and operators are taken at the scale /1 = Mw- 
At this point, it is convenient to use a Fierz transfor- 
mation. Such intrinsically four-dimensional identities are 
potentially problematic in d-dimensional renormalization 
schemes. In 't Hooft-Veltman dimensional regulariza- 
tion, the renormalization of operators is independent of 
whether they are Fierz transformed |34| . After a Fierz 
transformation, the charged-current contribution to the 
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isovector parity- violating Lagrangian has the form 



IV. TWO-LOOP RUNNING 



c 



'(5) 



Gj 



-IVcdfiuj^u-dj^Lic^L -{L^R} 



In writing the effective Lagrangian above 
dropped bottom quark contributions, 
cause they enter proportional to \V u i>Y 
2 x 10~ 4 is negligible. 



(29) 

have 
This is be- 
and the ratio 
Given the ex- 
, we can make 
0.2252(9), and 
This 



IK. 



\V ub \ 2 /\Vus? 

perimental values [39( for \V US \ and \V c d 
a further approximation. With \V US \ — 
\Vcd\ = 0.230(11), we shall take \V cd \ 
approximation is quite reasonable because the charged- 
current component of isovector hadronic parity violation 
is only a 10% contribution compared to the neutral- 
current component. Dropping the bottom quark oper- 
ators and treating \V cc i\ — \V US \ amount to neglecting 
sub-percent effects overall. 

Under these numerical approximations, we have the 
charged-current contribution to | A/| = 1 parity violation 
in the following form 



/(5) 



G F \V US \ 2 



V2 



djud] 



[ S7 M s - 



{L o R} 



(30) 



While the operator in Eq. (|3"0")) does not mix with those 



C 



(5) 



under QCD renormalization, an additional op- 



erator is required to form a closed set of charged-current 
contributions. The closed set has been given above in 
Eq. (J3J) - From comparing Eqs. (|3"0"j) and JSJ), moreover, 
we find the tree-level values of the remaining Wilson co- 
efficients in the five-flavor effective theory, 



A. QCD Renormalization Group 



Evolution of Wilson Coefficients 



The Wilson coefficients, c\ Ni \ijl) evolve with the scale 



H according to the renormalization group equation, 



L-j 3 



(a s ), 



(32) 



with T^ N f' as the anomalous dimension matrix in the 
theory with Nf flavors. The entry describes the 

renormalization of the operator O 
scribes its mixing into the operator O 



(N f ) 



while T^ j) de- 

J 1 



(N f ) 



It is conve- 



nient to drop the superscript (Nf) for notational ease. 
When we match between different effective theories, the 
superscript will reappear. 

At next-to-leading order, the anomalous dimension 
matrix can be written as 



71, 



(33) 



with the dependence on a s made explicit. Here 70 is 
the leading-order anomalous dimension matrix, and 71 
is the next-to-leading order matrix. The QCD running 
coupling satisfies the renormalization group equation 



-0o 



4-7T 



fa 



(4tt) 2 



(34) 



with fa = |(1L/V C - 2N f ), and fa = f A 2 - fN c N f 
(N c — N ( T 1 )Nf, the leading and next-to-leading order co- 
efficients of the /^-function, respectively. The solution 
to the renormalization group equation for a s (/j,) can be 
written as 



Us(lA 



-in 



/Jolog^/A 2 ) 



1 



fa loglog( M 2 /A 2 ) 
HI log( M 2 /A2) 



(35) 



which serves to define the QCD scale A. 

The solution to the renormalization group equation for 
the Wilson coefficients, Eq. (f32|) . can be written as 



cf\M w ) =0 
cf ) (M w ) = \V us ] C 



-0.0508. 



(31) 



Finally as the charged-current interaction leads to a 
subdominant contribution to |A7| = 1 hadronic parity 
violation, we do not consider its one-loop QCD correc- 
tions. Omitting such corrections amounts to neglecting 
sub-percent effects, which is consistent with the approx- 
imations used thus far. The tree-level values for Wilson 
coefficients will accordingly be used for charged-current 
contributions. 



Ci(n') 



(36) 



Here [/(//,//) is the evolution matrix. Valid to next-to- 
leading order, the explicit form of the evolution matrix 
is specified by [37J 



1 



47T 



J )U ( l i , ,fx) 1 



a s (fj,) 



4tt 



with Uq as the leading-order evolution matrix 



U (fj,',n) = V 



a 8 (n) 



(tq)-d 



V 



(37) 



(38) 
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defined in terms of the matrix V that diagonalizes 7^, 
namely (70)1? = V _1 7^V. The matrix J appears in the 
form 



where the matrix K is defined by 



2/3o + - [(7d)©]j 



(39) 



(40) 



2. Matching at Heavy Quark Thresholds 

The evolution of Wilson coefficients expressed in 
Eq. Q36[) is applicable as we run down from, say, fi = Mz 
to fi = m b . In order to evolve to lower scales, we must 
confront the matching from a five-flavor to a four-flavor 
effective theory, or, more generally, from an TV^-flavor to 
an (Nf — l)-flavor effective theory. Matching of theories 
at the heavy quark threshold has two parts: continuity of 
the coupling constant, and determination of the Wilson 
coefficients, c\ Nf (mo). 

As the QCD /3-function is Nf dependent, the cou- 
pling, a s (ji), will be discontinuous as we pass through 
a quark threshold, fi = vtiq. This unphysical discontinu- 
ity is avoided by adjusting the QCD scale parameter A 
in Eq. (|35l) . For a theory with Nf flavors, the QCD scale 
is then A^- Nf \ with continuity in the coupling enforced at 
the quark threshold, tuq, by fixing A^ -1 ) so that 

^''(mg.AW) = a£ w '- x W A^'- 1 *). (41) 

Having fixed the QCD coupling to be continuous, 
we must determine values for the Wilson coefficients, 
c\ Nf ^(itiq), from the c\ Nf \m.Q). In order to match 
the Nf- and (Nf — l)-flavored theories at the thresh- 
old fi = m,Q, we must first regroup the set of operators 

{o\ Nf ^} to match those in the (Nf — l)-flavor theory 
when the heavy quark is integrated out. It is best to 
spell this out explicitly for matching at the two scales 
/i = mb, and \i = m c . 

To match at the bottom quark threshold, we note that 
for Nf = 5, there are eight operators in the effective 
theory, and only six operators in the four-flavor theory. 

The third set of operators, 0y 5 g(/i), in Eq. ((3J) is already 

present in the four- flavor theory, namely as Orfl(p) in 
Eq. ©. As we integrate out the bottom quark, the val- 
ues of the Cy 5 g(m;,) coefficients will smoothly go over as 

contributions to the Cg 4 6 (mt). The remaining six opera- 
tors require matching conditions to avoid discontinuities 
when the bottom quark is integrated out. 

We now spell this matching out explicitly. At the 
bottom quark threshold, fi — m b , we integrate out 
the bottom quark at one-loop order in the five-flavor 
effective theory. This produces four-flavor operators: 



o[%(m b ), Of>(m b ) -> Omm b ), and 



,(5) 



7 g(m(,) = G§\(m b ). After regularization and subtrac- 
tion, the one- loop amplitudes are described by 0(a s ) fi- 
nite terms encoded in the matrix r^ 5 \ In a theory with 
Nf flavors, the matrix r"'' takes into account operator 
mixing at next-to-leading order, similar to r in Eq. (|19p . 
Here, the matrix r^ Nf ^ arises from the finite contribu- 
tions to the one-loop renormalization of all four-quark 
operators. 

Because the set of operators C 7 5 g (//) in Eq. ^ is closed 
under renormalization, the full matrix r^ 5 -* must have a 
block diagonal structure. The first block is an 6 x 6 di- 
mensional matrix, while the second block is a 2 x 2 di- 



mensional matrix. The latter matrix is identical to 



J4) 



for i, j = 5 or 6, because the corresponding operators are 
identical, and there is no effect at one-loop from integrat- 
ing out the bottom quark. In the five-flavor theory, the 
effective action at /i = m b is thus proportional to 



6 

E 



a s (rn b ) ( 5) 
4tt Tij 



0f>(m b ) 



2 

E 

i,j=l 



C\%{m b ) 



5 a + 



a s (m b ) ( 4 ) 



4tt 



i+4 J+4 



u 3 +i 



(42) 

while in the four-flavor theory, the next-to-leading order 
effective action at fi = m b is proportional to 



6 

E 



c\ A) (m b ) 



a s (m b ) ( 4 ) 
4tt 



of\m b ), (43) 



with the same constant of proportionality. In order that 
the theories match, we must have 



6 

E^: 

J'=l 



(5) 



(m b ) 



a s (m b ) (S ) 
4tt ji 

k&)C\%(m b ), 



+(5 i5 +S ie )C^ 2 (m b ), (44) 

where, in general, we have defined 

Ar^/) ee 7W - r^f-V. (45) 

Notice that to perform the matching at threshold, it is 
only the difference of the flavor-dependent mixing matri- 
ces, Eq. (|45|) . that is required. This difference incorpo- 
rates the finite one-loop contributions from integrating 
out the heavy quark. 

Matching at the charm quark threshold is simpler. At 
the scale fj, = m c , we integrate out the charm quark, and 
the operators o[ 4 \m c ) become o| 3 '(m c ), for each i = 
1, . . . , 6. This fact enables us to determine the matching 
condition, 



cfV c ) = £cf(m c ) 



a s (m c ) (4 ) 
4^ ^ 



(46) 



with the matrix Ar*- 4 ^ as the difference of finite contri- 
butions in the two effective theories, given in Eq. flU]) . 
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3. Solution: Wilson Coefficients at Hadronic Scales 

To summarize this section, we reiterate how to deter- 
mine the Wilson coefficients at hadronic scales, Aqcd < 

H < m c . The input values are C[_ 6 (Mz), and dfl(Mw) 
determined above in Sec. IIIII These values are evolved 
down to a common scale /i — mb, using the next-to- 
leading order evolution in Eq. (1551) . Fortunately these 
two sets of input values can be run down independently 
because they do not mix under renormalization. From 
the eight values C\ b ^ {mb), we match at the bottom quark 
threshold, via Eq. (|44[) . to determine the six Wilson 
coefficients in the four-flavor effective theory, C- (m&). 
These we run down to the charm quark mass, 



C?\m c ) 



3 



(47) 



Matching at /i = m c is accomplished using Eq. (|46j) 

(3) 

to determine C\ (m c ). The desired Wilson coefficients, 

/■q\ 

C\ (/x) , for Aqcd < M < m c, are then simply 

cfV) 



Y,U%\^m c )cf\m c ). (48) 

3 



Finally to perform the evolution from the weak scale 
down to hadronic scales, we need to know the anoma- 
lous dimension matrices T^ N f> to perform the running, 
and the difference of mixing matrices AA Nf ) to perform 
the matching. These quantities are discussed next. 



('0i7/iV'3)i(^/'37 Ml / ; 2)L- It will be helpful to keep the two 
classes of diagrams and two types of penguins in mind 
when discussing the renormalization. A thorough discus- 
sion of renormalization of four-quark operators at next- 
to- leading order is given in [30J. 



1. N f = 5 

Let us begin with the five-flavor effective theory valid 
at scales above the bottom quark mass. The operators 
Of ] are given in Eqs. ©, ©, and ©. In this theory, the 
anomalous dimension matrix, is best considered in 
blocks corresponding to these three sets of operators. As 
current-current renormalization can alter the color struc- 
ture of operators, but not the flavor structure, each block 
is closed under renormalization from current-current type 
gluon radiation. Due to their flavor structure, operators 
in the last block are not renormalized by QCD penguin 
diagrams. Thus these operators do not mix with the 
other two blocks. Operators in the second block can mix 
with the first block through penguin diagrams of type 
one. Finally operators in the first block are renormalized 
by both type one and two penguin diagrams, however, as 
these are the penguin operators themselves, their flavor 
structure is not altered outside their own block. Based 
on these general considerations, the anomalous dimen- 
sion matrix must have the block form 



p(5) = 



5P + 2P' 







P2x4 


^2x2 










C2x2 / 



(49) 



B. Anomalous Dimensions 

The operators, o[ Nf \ are renormalized by gluon radia- 
tion. There are two generic classes of Feynman diagrams 
according to which the gluon radiation can be classified. 
These diagram classes are current-current interactions, 
and QCD penguins. One-loop diagrams are depicted in 
Fig. [2J At two-loop order, one attaches an additional 
gluon in all possible places, and the diagrams can again 
be classified in these two categories. 

Current-current renormalization does not alter the fla- 
vor structure of the operator in question. QCD pen- 
guins, however, can lead to the generation of so-called 
penguin operators that contain the flavor-singlet com- 
bination proportional to *^2 q (jiq)v- At next-to-leading 
order, QCD radiation can also generate operators that 
contain ^2 q (qq)A- There are two species of penguins, 
moreover, depending on how the quarks in the opera- 
tor are contracted. Type-one penguins are formed from 
a closed quark loop, i.e. when spin contracted quarks, 
such as t/> 3 and ip 3 in O = {^il^^hi^^^L, are 
Wick contracted along with gluon interactions. Type- 
two penguins, on the other hand, are formed when 
quarks with uncoupled spin are Wick contracted along 
with gluon interactions, such as ^3 and ip 3 in O' = 



where the 4x4 blocks C, P, and P' are the anomalous 
dimension matrices from current-current renormalization 
and penguins of type one and two, respectively. The 
matrix C2X2, which encodes the anomalous dimensions of 
current-current radiation of Vl ® Vl — Vr ® Vr operators, 
is the upper-left 2x2 submatrix of C. The matrix P2x4 
consists of the top 2x4 submatrix of P and encodes 
the type-one pengiun radiation from Vl®Vl~ Vr <g> Vr 
operators. Each of these matrices has an expansion in 
a s , of which we consider only the first two terms, 



C 



4^ + fc) Cl < 

2 



p' - — p' + (—Y p' 



47T 



(50) 



so that 7q consists solely of the matrices Co, Pq, and 

Pq without any powers of a s , and ^y[ Nf ^ similarly consists 
only of C\, Pi, and P[. Furthermore as these matrices 
appear in the other effective theories, we postpone giv- 
ing the matrix elements until we have first discussed the 
general form in all three theories. The numerical fac- 
tors in front of penguins in Eq. (|49[) obviously arise from 
the number of flavor contractions that can produce such 
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contributions. In writing Eq. (1491) . we have made consid- 
erable use of the flavor-blindness of gluon interactions. 



2. N f = 4 

In the four-flavor effective theory, the operators are 
enumerated in Eqs. ([7} and ©. Treating these sets in 
blocks, each block undergoes current-current renormal- 
ization, but only the first block of operators can mix with 
QCD penguin operators. The anomalous dimension ma- 
trix in the four-flavor theory thus has the block-diagonal 
form 



r (4) 



c 



AP- 




2P' 

^2x2 



(51) 



where the matrices C, P, and P' are the same ones that 
appear in Eq. (|49|) . 



3. N f = 3 

Finally at hadronic scales, we have the three-flavor ef- 
fective theory. The relevant operators have been given 
in Eqs. (fit))) and (TTTT) . It is similarly efficacious to or- 
ganize the anomalous dimension matrix in blocks. Each 
block undergoes current-current renormalization as well 
as renormalization from penguin diagrams. As above, 
only the first block receives contributions from penguins 
of type two. In terms of the matrices C, P, and P' em- 
ployed above, the form of the anomalous dimension ma- 
trix must be 



r (3) 



C 



3P + 2P' 

Pz x 4 C2 x 2 



(52) 



the only relevant diagrams because they contain a heavy 
quark loop. 

In matching the reduced basis (i = 1, . . . , 6) of five- 
flavor operators at the scale \x = m\, to the four-flavor 
operators via Eq. (|4"4")l , we have the 6x6 dimensional 
matrix 



(54) 



The 4x4 block matrix rp contains the finite contributions 
from the one-loop mixing due to the type-one penguin 
contraction of the bottom quark, and rp 2x4 is the upper 
2x4 submatrix of rp. In matching the four-flavor theory 
to the three-flavor theory at the charm quark threshold 
via Eq. (146|) . we require the matrix 



Ar<*> = 



r P 



(55) 



which contains same type of contributions; this time due 
to the charm quark loop. Thus to perform the matching 
at heavy quark thresholds, it is sufficient to know the 
matrix rp. 



D. Matrix Elements 

We have seen that to specify the anomalous dimension 
matrices for each of the effective field theories, we must 
determine the three matrices C, P and P'. We do this to 
next-to-leading order. Furthermore, to match the effec- 
tive theories at heavy quark thresholds, we additionally 
need the matrix rp that arises from finite contributions 
to the heavy quark loops. 



C. Finite Mixing Matrices 

At next-to- leading order, one also requires the matrices 
r ( N f ) that characterize finite contributions to operator 
mixing at one-loop order. In each of the three effective 
theories, Nf — 5, 4, and 3, the general structure of the 
matrix is obviously identical to that of the rW 

deduced above. For example, in the four-flavor effective 
theory, we have 

r <«> = (-+V 2 -,,L). m 

where the subscripts denote the type of gluon radiation 
that leads to these finite contributions. As these r( N f' 
matrices are required to perform the matching of Wil- 
son coefficients at heavy quark thresholds, moreover, we 
only need their differences, Ar^ Nf ^ defined in Eq. (|45|). 
At one-loop order, the only dependence on the number 
of flavors is generated by the type-one penguin contrac- 
tion, see Eqs. (l4*9")l . (|5"T1) . and (|52"1) . Physically these are 



1. Current- Current Anomalous Dimensions 

The matrix C arises from renormalization of the 
current-current interaction. We use a mass-independent 
renormalization scheme, so that the quarks are massless 
in each effective theory. Additionally the QCD radiation 
is flavor blind. These two features allow us to determine 
C from the renormalization of operators with four distin- 
guishable quarks, ipi, ...,ip±. In order to correspond to 
our 4x4 blocks, these operators are written as 

AQi = (?i7^2)l(^ 3 7^ 4 )l - {L O R}, 
AQ 2 = (?i7m^ 2 ]l[^ 3 7 m V'4)l - {L o R}, 
AQ 3 = (?i 7^2)^37^4)*. - {L R}, 
AQ 4 = (^^Wsl^^R -{L^ R}. (56) 

Denoting the anomalous dimension matrix of these oper- 
ators as Taq, we have simply C = T&q. 

Next we observe that parity transformed operators 
share the same anomalous dimension. This fact owes to 
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the parity invariance of QCD, as well as that of the reg- 
ularization scheme. Finally in the augmented 't Hooft 
Veltman scheme at next-to-leading order, the Vl <& Vl 
and Vl ® Vr operators do not mix under renormalization. 
Their anomalous dimension matrices, Tq vll and ro VLR , 
can be identified from the results determined in [34| . In 
that work, care was taken to handle evanescent opera- 
tors that appear at intermediate stages of the calculation. 
From these matrices, we form 



o r Q 



o 

VLR 



and consequently find 



^2x2 — Tc 



(57) 



(58) 



For completeness, the required anomalous dimension 
matrices for Vl ® Vl operators are 



(7Qvll)o 



-2 6 N 
6 -21 ' 



at leading order, and 

(7Qvll)i : 



Sf-fN, f-2N f 



2 ZiV / 6 9 iV // 



(59) 



(60) 



at next-to-leading order [34] . As for the Vl ® Vr opera- 
tors, we have the anomalous dimension matrix 

(7Qvlr)o = fg fgj , (61) 



at leading order, and the matrix 



121 -&N t -39-4N 



(7Qvlr)i = 

at next-to-leading order [3 



3 iV / 

2 3 J 2 9 Iy J/ 



(62) 



2. Penguins of Type One 

Renormalization of QCD penguins of type one leads 
to the anomalous dimension matrix P employed above. 
This matrix can be determined by studying how four 
generic operators are renormalized by QCD penguins of 
the first type. Due to flavor blindness of the QCD inter- 
action, these operators can be chosen as 



AP t = AQ 4 



(63) 



with AQi shown in Eq. ([51)1) . Notice that by design, the 
APi operators only have penguin contractions of type 
one. 

The renormalization of the chiral-basis operators, APi, 
from penguin diagrams can also be determined from the 



results of [3J], owing to the parity invariance of QCD and 
of the regularization scheme. To determine the matrix 
P, it is sufficient to study how the set of four operators 
{APi}, mixes with the penguin operators, 

N f 

APi = (^7^2)1, ^(97 m 9)l - {L O R}, 

9 
Nf 

AV 2 = (?i7„V>2]l Y}^^ l ~ i L fl >< 

9 

N f 

AV 3 = (Ai^2)l J2^^x ~{L^ R}, 

Q 

N f 

AV A = (Vi7^ 2 ]l Yfi^fiR ~ < L ° R y ( 64 ) 
i 

This mixing is described by a 4 x 4 matrix Tap, and our 
desired matrix is simply P = Tap- 

For completeness, the leading-order type-one penguin 
anomalous dimension matrix has the form 



(7ap)o 



/ 0\ 

2 2 _2 2 

"9 3 9 3 



2 2 _2 2 

" 9 3 9 3 



(65) 



While at next-to-leading order, the required anomalous 
dimension matrix can be deduced from the results of pjij . 
We find 

23 1 25 



/ f 1 -f 1\ 



(7ap)i = 



418 
" 243 



73 
9 



850 1210 490 
81 243 81 



71 

9 



(66) 



, 1246 382 256 832 / 
' 243 81 243 81 ' 



Considerable care was taken in [3J| to ensure that con- 
tributions from evanescent and off-shell operators were 
taken into account at intermediates stages of the calcu- 
lation. 

The final type-one penguin contributions we require 
are their finite contributions at one-loop. These contri- 
butions are contained in the matrix 37] 



/o o\ 

_5_ _5 _5_ _5 

27 9 27 9 



I _5_ _5 _5_ _5 i 

\27 9 27 9/ 



(67) 



As deduced above, this is the only matrix of finite con- 
tributions necessary to match effective theories at heavy 
quark thresholds. 



3. Penguins of Type Two 

To deduce the form of the anomalous dimension ma- 
trix that encodes the mixing of operators due to QCD 
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penguins of type two, we appeal to flavor blindness to 
consider the set of four operators, 



AP' t = AQi 



'/ , 2=V , 3,V'4=V'2 



(68) 



that, by design, has only the requisite penguin contrac- 
tions of type two. The mixing of AP( operators with the 
penguin operators, AVi, in Eq. (|64|) is described by an 
anomalous dimension matrix r^p', which is our desired 
matrix, P' = Tap' • 

The use of Fierz identities is valuable here. In 't 
Hooft-Veltman dimensional regularization, the two-loop 
anomalous dimensions are the same between Fierz trans- 
formed operators. The operators in Eq. (|68]) are related 
by a Fierz transformation to operators with only type-one 
penguins. In essence, we are using the reverse of the ar- 
gument employed by [37j to determine penguin contribu- 
tions. The Fierz transformation of the Vl <S> Vl — Vr <8> Vr 
operators produces two operators that we explicitly con- 
sidered above in dealing with type-one penguins, namely 



AP{ = AP 2 , 
AP2 = AP X . 



(69) 



On the other hand, the Fierz transformation of the two 
Vl®Vr— Vr®Vl operators are themselves new operators, 



AP; = 2 0M 2 ]l [fafoU - {L O R}, 
API = 2 QMl^s^r -{L^ R}, 



(70) 



that involve only scalar bilinears. Because we utilize a 
mass-independent renormalization scheme, closed quark 
loops arising from scalar bilinears will always involve an 
odd number of Dirac matrices, and thereby vanish. Thus 
the lower two rows of P' must be zero. 

From the results quoted above for AP 12 , we can eas- 
ily determine the type-two penguin anomalous dimension 
matrix. For completeness, we display these matrices ex- 
plicitly. The leading-order matrix has the form 



(TAP' 



2 


2 


2 


2 


9 


:i 


9 


3 




































(71) 



While at next-to-leading order, we have the matrix 



(7AP')l = 



/- 


418 


.sr.u 


1210 


490 \ 
81 1 


243 


81 


243 




23 


1 


25 


1 




3 


3 
















V 











/ 



(72) 



4- Consistency Checks 

As a consistency check, one can easily verify that after 
a linear transformation, the matrices Co, Po, and Pq lead 



to the same one- loop anomalous dimension matrix Jq N 
employed by [29| . 

A further consistency check is provided by the redun- 

(2) 

dancy of the operator 0\ expressed by the Fierz relation 
in Eq. (|12[) . To expose this redundancy, we can transform 



to a new basis specified by O 



'(3) 



Of, with the 



transformation mati 


ix 


jiven by 


























1 


-1 








-1 


1 


S = 








1 




















1 






















1 







Vo 














1/ 



(73) 



Thus the primed basis is identical to the unprimed basis 

(3) 

except for the redundant operator 0\ ' . In the primed 
basis, the corresponding operator is 



O 



'(3) 



(3) 



o 



(3) 



o 



(3) 



o, 



(3) 



(74) 



— . / (3) 

The Fierz identity implies 2 = 0, and accordingly 

this operator should not mix with any other operators. 

Remember that the 't Hooft-Veltman scheme preserves 

Fierz relations. 

In terms of the anomalous dimension matrix in the 

primed basis, r'( 3 \ we must have the vanishing of the 

0, for any i ^ 2. This condition 



components T^ 3 ' 



ensures that whatever the value of C 2 (//), it will have 
no effect on the evolution of the other Wilson coefficients, 
C^ 3 2(M),seeEq. fl32D0 The anomalous dimension matrix 

in the primed basis is given by r'( 3 ' = ST^S~\ where 
it happens that = S. We have the constraint 



ST (3) S 



2 j 



0, for i + 2, 



(75) 



which must be satisfied order-by-order in perturbation 
theory, and the leading-order and next-to-leading order 
anomalous dimension matrices indeed satisfy Eq. (|75l) . 

To next-to-leading order, we actually find that current- 
current contributions to Eq. (|75p vanish independent of 
the value of Nt (recall that C\ depends on Nf). The 
remaining contributions to the constraint involve the 
anomalous dimension matrices P and P' with particular 
numerical coefficients depending on Nf. For these terms, 
only the contributions for the value Nf — 3 vanish, and 
this is precisely the situation where the Fierz constraint 
exists. 



6 While the Wilson coefficient evolves with the renormal- 



-C. 



'(3) 



^'(3) r '(3) 



^ 0, the 



ization scale according to fj,-^- 
product C' 2 (S) (u)0' 2 "" 

the Fierz identity implies Q' 2 ^°' (/j,) = 0, at all scales. Conse- 
quently any matrix elements of this operator computed by lattice 
techniques should vanish up to statistical noise. 



' ^ (u)C>2^ (ft) is trivially scale independent because 
V(3)/ 
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E. Results 

With the anomalous dimension matrices fully specified, 
we can now run the Wilson coefficients down to hadronic 
scales. We carry out the running in the chiral basis; and, 
for the reasons described above in Sec. IIII1 we find is 
practical to convert to the |A/| = 1 parity-violating op- 
erator basis given in [l8[ . This basis is specified by eight 
operators. There are four non-strange operators 

0\ = (uu — dd)A(uu + dd)y, 

02 = (uu — dd]A[uu + dd)y, 

03 = (uu — dd)v(uu + dd)^, 

04 = (uu — dd ] v [ uu + dd)A, (76) 

and four strange operators 

05 = (uu — dd)A(ss)v, 

6 = (uu - dd]A[ss)v, 
Oj = (uu — dd)v(ss)A, 

8 = (uu — dd] v [ss)A- (77) 

To convert between bases, it is convenient to append 
two operators to Eq. (fTT|) . namely 

Of ] = (u 7/i u - d^d) L (s^s) R -{L^R} 

( s 3) = (u 7 „u - dlnd} L [s"(> l s) R -{L^ R}. (78) 

These operators necessarily have zero Wilson coefficients. 

The set of operators Oi is related to the set of eight op- 

(3) 

erators 0\ through a series of linear transformations, 

jk 

with the transformation matrix T responsible for con- 
verting the chiral basis to the vector-axial- vector basis. 
It is specified by T = diag(T, T) , with the 4x4 transfor- 
mation matrix T given by 

T =~(] -J)' ( 8 °) 

in 2 x 2 block form. The transformation matrix J alters 
the flavor structure of the operators so that we arrive at 
a separation between strange and non-strange operators. 
Written as a 4 x 4 block matrix, we have 

■ 7 =(J"1)- (81) 

Finally, for the eight Wilson coefficients in the new basis, 
we write C;(/i), and these are given by 

£ t (u)=J2c { k 3) (v)T kj 1 J k ~\ (82) 

jk 

with c4 3 8 V) = 0. 



TABLE I. Values for the Wilson coefficients of the |A/| = 1 
parity- violating operators, £i(u), at hadronic scales. We 
compare values obtained from leading-order (LO) evolution 
quoted in to results obtained in this work. In all cases, 
we normalize the coefficients by dividing by the tree-level co- 
efficient j(co)i| = I sin 2 8w, and the values quoted are at the 
scale jj, = 1 GeV. Ingredients of our LO and next-to-leading 
order (NLO) computations are described in the text. To make 
a meaningful comparison of the results, we eliminate the re- 
dundant operator O4 using Eq. (183(1 . 







Ci(lGeV)/|(co)i| 




i 


LO [18] 


LO 


NLO (Z) 


NLO (Z+W) 


1 


0.403 


0.264 


-0.054 


-0.055 


2 


0.765 


0.981 


0.803 


0.810 


3 


-0.463 


-0.592 


-0.629 


-0.627 


4 














5 


5.61 


5.97 


4.85 


5.09 


6 


-1.90 


-2.30 


-2.14 


-2.55 


7 


4.74 


5.12 


4.27 


4.51 


8 


-2.67 


-3.29 


-2.94 


-3.36 



Notice that in the basis employed by [l8|], the Fierz 
transformation leads to the constraint 

4 = O x - 2 + 3 . (83) 

Finally as a consequence of the renormalization in the 
chiral basis, the operators 

0\ = -O x + 3 + 5 - 7 

0' 2 = -0 2 + O i + 6 -O s , (84) 

have vanishing Wilson coefficients because they are not 
generated by QCD radiative corrections. 

Computed values for the Wilson coefficients £i (/i) at a 
scale of /i = 1 GeV are collected in Table HI As inputs to 
this calculation, we use the masses |39j |. 

M z = 91.2 GeV, M w = 80.4 GeV, 

m b = 4.19GeV, m c = 1.29GeV, (85) 

where the latter two are MS masses, and correspond to 
the masses in the subtraction scheme we employ. The 
value of the strong coupling is taken at the weak scale, 

a s (M z ) = 0.118. (86) 

From the two-loop expression for a s ([i) given in Eq. (|35l) , 
we find the value of the five-flavor QCD scale parameter 
A^ 5 ) = 0.231 GeV. The requirement of continuity at nib 
and m c leads to the values A( 4 ) = 0.330 GeV and A^ 3 ) = 
0.376 GeV, respectively. 

The results shown in the table compare our leading and 
next-to-leading order computations of the renormaliza- 
tion of isovector parity-violating operators. Ingredients 
of these computation are summarized as follows: 
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1. Leading order (LO): coefficients evolved with one- 
loop anomalous dimensions but with two-loop run- 
ning coupling, tree-level matching at the weak scale 
and at heavy quark thresholds, charged-current in- 
teractions excluded. 

2. Next-to-leading order Z [NLO (Z)]: coefficients 
evolved with two-loop anomalous dimensions along 
with the two-loop running coupling, next-to- 
leading order matching at the weak scale and at 
heavy quark thresholds, charged-current interac- 
tions excluded. 

3. Next-to-leading order Z+W [NLO (Z+W)]: same 
as NLO (Z), but now including charged-current in- 
teractions at tree level. 

We also compare our values with the lead ing -order evolu- 
tion of Wilson coefficients determined in [29j , and quoted 
at /i = 1 GeV in [l8[. To present these and our results, 
moreover, we do not give values for all eight coefficients. 
The Fierz constraint on the operators, Eq. (|53")) . makes 
reporting values for the first four coefficients somewhat 
arbitrary. To eliminate this ambiguity, we choose O4 as 
the redundant operator, and accordingly determine con- 
tributions to £i-3([i) from its coefficient £a{^)- 

In the table, we see that from comparing leading and 
next-to-leading order results, the corrections are gener- 
ally on the order of 10-20%. This is consistent with the 
typical size of next-to-leading order effects in the context 
of the QCD evolution of flavor-changing weak currents. 
We also see that inclusion of isovector charged-current 
interactions has little effect on the coefficients of non- 
strange operators: the first three coefficients are altered 
only by a few percent when the charged-current contri- 
butions are turned on. This behavior is not surprising 
given that the charged-current operators only mix with 
non-strange operators through penguin diagrams, and 
this happens in the QCD evolution, moreover, only when 
fi < m c . Were we to evolve the coefficients down to lower 
scales, charged-current interactions would have a larger 
effect on the non-strange operator coefficients. The ef- 
fect of charged-current interactions on the coefficients of 
strange operators, by contrast, is at the 10-20% level. 

Finally, the behavior of the coefficient gives one 

reason to pause. Our value at leading order is quite dif- 
ferent from that in [Lsj , which is already an indication of 
its sensitivity to higher-order corrections. The difference 
between leading and next-to-leading order is 120% of the 
leading-order value. While the size of such corrections 
is alarming, this coefficient is the smallest (in absolute 
value) of all hadronic-scale coefficients. This is true even 
if we use the leading-order value for Assuming the 

isovector parity- violating hadronic matrix elements of all 
Oi (/it) operators are roughly the same size, we should not 
compare the difference of leading and next-to-leading or- 
der contributions to the value of GTi (/x) , per se, but rather 
to the maximum value of all the coefficients. The 

next-to-leading order corrections to are only 6% 



of the maximum Wilson coefficient. Overall the next- 
to-leading order corrections to isovector parity-violating 
hadronic couplings should be under control. This can be 
spoiled in practice if parity- violating hadronic matrix el- 
ements of Oi(fi) are dynamically enhanced an order of 
magnitude over the other operators, or if there are dra- 
matic cancelations among the other terms. 



V. OUTLOOK 

Above, we study the QCD evolution of isovector 
parity- violating four-quark operators from the weak scale 
down to hadronic scales. At hadronic scales, matrix 
elements of the corresponding operators are required 
to determine the |A/| = 1 parity-violating couplings 
between hadrons non-perturbatively using lattice QCD 
techniques. We find that next-to-leading order contribu- 
tions alter the values of Wilson coefficients at fi = 1 GeV 
by ~ 10—20% in most cases. To compute these values, we 
determine the one- loop matching of Wilson coefficients at 
the weak scale, Eq. (121)1) , and deduce the form of isovec- 
tor parity violation from the charged-current interaction, 
Eq. pip . Appealing to parity invariance and the flavor- 
blindness of QCD radiation, we extract the required two- 
loop anomalous dimension matrices from those arising in 
the study of QCD corrections to |A5*| ~ 1 non-lcptonic 
weak decays |34| . Results are given in the dimensional 
regularization scheme of 't Hooft-Veltman, with an addi- 
tional multiplicative renormalization of the axial- vector 
current which is necessary in order to preserve chirality. 
In this scheme, we find that next-to-leading order cor- 
rections to Wilson coefficients are important. Moreover, 
inclusion of |A/| = 1 charged-current interactions affects 
the coefficients of strange operators at hadronic scales by 
10-20%. 

Having determined the Wilson coefficients of isovec- 
tor parity- violating operators at hadronic scales to next- 
to-leading order accuracy, there are a number of fur- 
ther studies that are now open to pursuit. First is the 
inclusion of electromagnetic corrections. In this work, 
we focus solely on QCD corrections; and, at hadronic 
scales, fi, we include corrections that are of the order 
[a s (/x)/(47r)] 2 . The leading electromagnetic corrections 
arise proportional to a/4ir, and are thus expected to be 
of the same size. Further study is required to determine 
the effects of electromagnetism on parity-violating op- 
erators at hadronic scales, and this involves mixing be- 
tween the various isospin channels. A second investiga- 
tion is the study QCD renormalization of parity viola- 
tion using a lattice regularization. One-loop computa- 
tions using lattice perturbation theory are required to 
determine the renormalization scheme dependence, and 
knowledge of this result will enable one to convert a 
lattice regularization to the MS scheme. Such knowl- 
edge is essential to make physical predictions from lat- 
tice QCD computations of hadronic parity violation. To 
this end, a final investigation is to study the evolution of 
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Wilson coefficients using a non-perturbative renormaliza- 
tion scheme, such as the renormalization independent — 
momentum scheme [loi. Nevertheless, the investigation 
of hadronic parity violation represents an exciting op- 
portunity for lattice QCD. In light of major experimen- 
tal efforts, the lattice approach will be complimentary 
in mapping out the parity- violating nuclear force, and 
will ultimately connect these few-body parameters to the 



hadronic weak interaction in the Standard Model. 



ACKNOWLEDGMENTS 

Work supported in part by a joint CCNY-RBRC fel- 
lowship, and by the Research Foundation of the CUNY. 
We thank J. Wasem for enlightening discussions, and the 
hospitality of Lawrence Livermore and Lawrence Berke- 
ley National Laboratories during visits in the course of 
this work. 



[1] S. Glashow, J. Iliopoulos, and L. Maiani, Phys.Rev. D2, 

1285 (1970) [22 

N. Tanner, |Phys.Rev. 107, 1203 (1957) 

V. Lobashov, V. Nazarenko, L. Saenko, L. Smotritsky, [23 
and G. Kharkevitch, |Phys.Lett.| B25, 1 04 (1967) 

C. Bowman, J. Bowman, and V. Yuan, Phys.Rev. C39, [24 
1721 (1989) 

E. Adelberger and W. Haxton, Ann. Rev. Nucl. Part. Sci. [25 
35, 501 (1985) 

M. J. Ramsey-Musol f and S. A. Page, [26 
I Ann. Rev. Nucl. Part. ScTl 56, 1 (2006), [27 
|arXiv:hep-ph/0601127 [hep-ph]] 

W. M. Snow et al., |Phys.Rev.|C 83, 022501 (2011) [28 
M. T. Gericke et al., |Phys.RevT] C83, 015505 (2011) 
B. Desplanques, J. F. Donoghue, and B. R. Holstein, [29 
|Annals PhysJ, 124, 449 (1980) 

S.-L. Zhu, C. Maekawa , B. Holst ein, M. Ramsey- [30 
Musolf, and U. van Kolck,|Nucl Phys.| A748, 435 (2005), 
|arXiv:nucl-th/0407087 [mid-th]] - 

D. R. Phil lips, M. R. Sch indler, and R. P. Springe r, [31 
|Nucl.Phys. A822, 1 (2009), |ajXTv:0812.2073 [nucl-th]| [32 
P. F. Bedaque and U. van Kolck, Ann. Rev. Nucl. Part. Sci. 
52, 339 (2002), |arXiv:nucl-th /0203055 [nucl-th] | [33 
L. Girlanda, |Phys.Rev.| C77; 067001 (2008), 

|arXiv:0804.0772 [nucl-thJT [34 

J. Shin S. Ando, and C Hyun, |Phys.Rev. C81, 055501 
(2010), |arXiv:0907.3995 [nucl-th]] 

M. R. Sch mdler and R. P. Springer , |Nucl.Phys.| A846, [35 
51 (2010), |arXiv:0907.5358 [nucl-th]| 
H. W. Griesshammer, M. R. Schindler, and R. P. 
Springer(2011), |arXiv:1109.5667 [nucTthJI [36 
J. Vanasse(2011), |arXiv:1110.1039 [nucbthj] [37 
D. B. Kaplan and M. J. Savage, |Nucl.Phys.| A556, 653 
(1993) [38 
S.-L. Zhu, S. Puglia, B. R. Holstein, and 
M. Ramsey-Musolf, |Phys.Rev.| D63, 033006 (2001), [39 

|arXiv:hep-ph/0005281 [hep-ph] | 

P. F. B edaque and M. J. Savage, |Phys.Rev. |C62, 018501 [40 

(2000) , |arXiv:nucl-th/9909055 [nucl-th] | 

J.-W. Chen, T. D. Cohen, and C.-W. Kao, |Phys.Rev.| 



C64, 055206 (2001), |arXiv:nucl-th /0009031 [nucl-th] 
J.-W. Chen and X.-D. Ji,|Phys.Rev.Lett. 86, 4239 (2001), 
|arXiv:hep-ph/0011230 [hep-ph] | 

J.-W. Chen and X.-D. Ji, |Phys.Lett.] B501, 209 (2001), 

|arXiv:nucl-th/0011100 [nucl-th] | 

S. R. Beane and M. J. Savage, |Nucl.Phys.| B636, 291 

(2002), |arXiv:hep-lat/020 3028 [he p-lat]| 

B. C. Tiburzi, |Phys.Rev.| D71, 054504 (2005), 



|arXiv:hep-lat/0412025 [hep-lat] | 

J. Wasem(2011), |arXiv:1108.1151 [heplat]] 
T. Blum, P . Boyle, N. Christ, N Garron, E. Goode, et al, 
|Phys.Rev . D84, 114503 (2011), [arXTv:1106.2714 [hep-lat]| 
T. Blum, P. Boyle, N. Christ, N. Garron, "E~Goode, 



Laut- 
(1996), 



et al. (201 1) , [arXiv: 1 1 1 1 . 1699 [hep-lat]] 
J. Dai, M. 3. Savage, T. Liu, and R. P. Springer, 
|Phys.Lett.| B271, 403 (1991) 
G. Buchalla, A. J. Buras, and M. E. 
enbacher, |Rev.Mod.Phys.| 68, 1125 
|arXiv:hep-ph/9512380 [hep-ph] | 

G. 't Hooft and M. Veltman, |Nucl.P"h7sTlB44 ,189 ( 1972) 
D. Akyeampong and R. Delbourgo, Nuovo Cim. A17, 
578 (1973) 

P. Breitenlohner and D. Maison, Commun.Math.Phys. 
52, 11 (1977) 

A. J. Buras, M. Jamin, M. E. Lautenbacher, 
and P. H. Weisz, |Nucl.Phys.|| B400, 37 (1993), 
|arXiv:hep-ph/921 1304 [hep-ph] | 

Ciuchini, E Franco, Lubicz, G. Martinelli, 
I. Scimemi, et al, |Nucl.Phys.| B523, 501 (1998), 
|arXiv:hep-ph/9711402 [hep-ph]| 

A. J. Buras and P. H. Weisz, |Nucl.Phys.| B333, 66 (1990) 
A. J. Buras, M. Jamin, M. Lautenbacher, and P. H. 
Weisz, Nucl.Phys.| B370, 69 (1992) 

M. Ciuchini, E. Franco, G. Martinelli, and L. Reina, 
|Nucl.Phys.| B415, 403 (1994), [hep^ph/9304257 
K. Nakamura et al. (Particle Data Group), 
G37, 075021 (2010) 
G. Martinelli, C. Pit tori, C. T. Sachrajda, M. Testa, 
and A. Vladikas, |Nucl.Phys.| B445, 81 (1995), 
|arXiv:hep-lat/9411010 [hep-lat] | 



J.Phys.G 



